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ANALYTICAL  INVESTIGATION  OF  HIGH  PERFORMANCE 
SHORT,  THRUST  AUGMENTING  EJECTORS 


Tah-teh  Yang  and  Francois  Ntone 
Mechanical  Engineering  Department,  Clemson  University 
Clemson,  South  Carolina 

ABSTRACT 

A  procedure  for  analyzing  thrust  augmenting  ejectors  using  a 
short,  curved-wall  diffuser  is  presented.  In  this  type  of  diffuser 

a  shear  flow  is  admitted  at  the  inlet  and  a  set  of  auxiliary 
ejectors  is  used  to  provide  the  necessary  boundary  layer  control. 

Several  computer  programs  are  used  in  the  analytical  investigation. 
Analyses  of  these  computer  programs  are  either  outlined  or 
referenced  in  the  open  literature.  A  user’s  manual  is  provided  in 

the  appendices  of  this  report.  A  discussion  of  the  geometries  and 

performances,  including  thrust  augmentation  ratios,  of  two  sample 
ejectors  determined  by  use  of  the  present  analysis  procedure  with 
mixing  chamber  contraction  as  an  optimization  parameter  is 
presented.  Use  of  the  mixing  chamber  length  as  an  optimization 
parameter  is  also  discussed. 


ADMINISTRATIVE  INFORMATION 

THE  ANALYTICAL  WORK  PRESENTED  IN  THIS  REPORT  WAS  SUPPORTED  BY  THE  DAVID 
TAYLOR  NAVAL  SHIP  RESEARCH  AND  DEVELOPMENT  CENTER  AND  NAVAL  AIR  SYSTEMS 
COMMAND  UNDER  CONTRACT  DTNSRDC  N00I67-80-C-0040,  TECHNICALLY  UNDER  THE  COGNI¬ 
ZANCE  OF  D.  KIRKPATRICK  (NAVAIR  320D)  AND  T.  C.  TAI  (DTNSRDC  1606). 


1 .  INTRODUCTION 


The  thrust  augmentation  ratio  of  an  ejector  depends  to  a  large  degree  on 
the  effectiveness  of  its  diffuser  (1,  2,  3).  In  early  1970  Yang  and  El-Nasher 


(4)  carried  out  an  exploratory  investigation  with  an  ejector  having  a  short 
curved-wall  diffuser.  They  used  saturated  steam  for  the  primary  flow  and 
ambient  air  as  the  secondary  flow.  Their  objective  was  to  examine  the  poten¬ 
tial  increase  in  ratio  of  the  mass  flow  rate  at  the  ejector  exit  to  the 
primary  flow  at  the  nozzle  exit.  The  pressure  ratio  range  of  their  experi¬ 
ments  was  40  to  400*.  On  the  average  the  mass  ratios  were  more  than  twice 
those  reported  by  Hickman  (5).  In  Yang  and  El-Nashar's  experiments  a 
centrifugal  pump  provided  the  boundary  layer  control  for  the  short  curved-wall 
diffuser,  and  the  rate  at  which  the  fluid  was  removed  from  the  boundary  layer 
was  less  than  15%  of  the  total  flow  rate  at  the  exit  of  the  ejector  or  the 
diffuser.  Based  on  the  mass  flow  rate  and  the  velocity  profiles  at  the 
diffuser  exit,  a  maximum  thrust  augmentation  ratio  of  2.81  was  estimated.  The 
high  value  of  the  thrust  augmentation  ratio  is  partly  attributed  to  the  use  of 
steam  as  the  primary  flow.  This  value  would  have  been  approximately  2  if  air 
of  the  same  stagnation  pressure  but  of  a  lower  stagnation  temperature  had  been 
used  for  the  primary  flow.  However,  the  most  important  feature  of  this 
ejector  is  its  overall  length.  Because  a  short  curved-wall  diffuser  was  used, 
the  length  of  the  ejector  was  reduced  to  about  one-half  the  length  of  a 
conventional  ejector.  This  reduction  in  length  is  particularly  attractive  for 
application  where  the  length  restriction  is  severe.  The  disadvantage  of  such 
an  ejector  is  the  need  for  boundary  layer  control.  A  scheme  for  using  a  set 
of  secondary  nozzles  to  form  an  auxiliary  ejector  and  thus  to  provide  the 
necessary  boundary  layer  control  was  evolved.  The  discharge  from  this 
auxiliary  ejector  is  oriented  in  the  same  direction  as  the  main  ejector  and 
therefore  contributes  to  the  overall  thrust;  this  makes  the  use  of  boundary 
layer  control  more  acceptable.  Figure  1  shows  a  sketch  of  a  thrust  augmentor 
equipped  with  a  short  curved-wall  diffuser  where  a  set  of  secondary  jets 
provide  boundary  layer  suction  for  the  diffuser. 

The  short  curved-wall  diffuser  was  developed  at  Clerason  University  during 
the  late  1960s  and  early  1970s.  Details  of  the  analytical  and  experimental 
flow  fields  within  these  diffusers  have  been  reported  in  engineering  journals 


*  Following  jet  pump  convention,  this  pressure  ratio  is  defined  as 
(P0l-Pe)^pe-pl>* 
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(6,  7,  8,  9,  10),  as  well  as  in  contractor  reports  (11,  12,  13,  14).  These 
diffusers,  designed  by  using  the  Inverse  Method,  have  sometimes  been  referred 
to  in  the  literature  as  Griffith  diffusers  because  the  velocity  distribution 
along  the  diffuser  wall  is  taken  from  the  high  lift  airfoil  designed  by  A.  A. 
Griffith  in  1938.  A  typical  Griffith  Diffuser  wall  velocity  distribution  is 
shown  in  Figure  2.  It  exhibits  a  high  velocity  inlet  region,  slightly 
accelerating  flow  up  to  a  narrow  region,  and  a  sudden  drop  in  velocity  across 
a  narrow  slot  region.  Downstream  of  the  slot  region  there  is  an  exit  region 
toward  Che  end  of  the  diffuser.  Figure  3  shows  the  typical  shape  of  the 
diffuser  which  is  the  result  of  the  prescribed  wall  velocity  distributions. 
Figure  4  shows  several  velocity  distributions  of  the  potential  flow  along  the 
streamlines  away  from  the  diffuser  wall. 

Since  there  is  no  deceleration  over  the  solid  surface,  flow  separation 
cannot  take  place.  However,  if  fluid  is  to  jump  over  the  slot  region,  it  must 
have  a  velocity  higher  than  the  critical  velocity.  The  critical  velocity  can 
be  computed  by  using  Taylor's  criterion  (15).  Much  experimental  data  have 
been  examined  for  diffuser  inlet  flows  having  thin  boundary  layers  with  core 
flows  well  represented  by  potential  flow  theory.  In  such  cases,  experimen¬ 
tally  obtained  wall  pressure  distributions  agree  with  the  theoretically  pre¬ 
scribed  distributions.  In  ejectors,  however,  the  diffuser  inlet  flow  is  sig¬ 
nificantly  different  from  irrotational  flow. 

Based  on  the  inlet  velocity  measurements  reported  by  Hill  and  Gilbert 
(16),  it  appears  that  the  diffuser  inlet  flow  should  be  represented  by  a  shear 
flow  of  uniform  vorticity.  Under  this  condition  (admitting  a  shear  flow)  the 
wall  pressure  distribution  could  be  significantly  different  from  that  which 
was  prescribed  for  the  particular  diffuser.  Specifically,  this  raises  a  con¬ 
cern  about  the  presence  of  an  adverse  pressure  gradient  or  a  deceleration  in 
diffuser  wall  velocity  which  may  result  from  the  inlet  shear  flow.  A  moderate 
deceleration  upstream  of  the  suction  slot  may  be  overcome  by  increasing  the 
suction  flow  rate.  A  strong  deceleration  will  result  in  flow  separation.  The 
combination  of  a  large  design  area  ratio  and  inlet  vorticity  could  also  create 
a  flow  reversal  at  the  diffuser  exit  because  of  the  shear  flow  requirement, 
even  in  an  inviscid  flow.  it  becomes  apparent,  therefore,  that  a  method  of 
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calculating  the  wall  velocity  of  a  rotational  flow  in  this  type  of  curved 
diffuser  is  necessary  to  assure  that  there  will  be  no  adverse  pressure 
gradient  along  the  diffuser  wall.  This  will  allow  us  to  obtain  the  high 
thrust  augmentation  thrust  augmentation  benefit  afforded  by  using  a  short, 
curved-wall  diffuser  section  in  the  ejector  design. 


2 .  ANALYSES 


The  analyses  in  this  investigation  will  be  presented  in  the  order  out¬ 
lined  in  Figure  5. 

2.1  MIXING  CHAMBER  GEOMETRY 


The  geometry  of  a  mixing  chamber  is  described  by  the  equation 

C 


R 


i  1  — 

in, me  \  R_. 


2  3 

10,x/' 


(1) 


in, me 

where  r  is  the  radius  at  any  location  of  x 


Itin.mc  is  the  radius  of  the  mixing  chamber  at  the  inlet 
x  *  0  is  the  location  of  the  mixing  chamber  inlet 
x  *  L  is  the  location  of  the  mixing  chamber  exit 

C  is  a  dimensional  constant. 

Thus,  contraction  is  gradual  and  the  maximum  contraction  takes  place  at 
the  exit  of  the  mixing  chamber.  This  geometry  provides  zero  wall  slopes  at 
the  inlet  and  exit  of  the  mixing  chamber. 
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2.2  DTNSRDC  EJECTOR  PROGRAM 


The  DTNSRDC  ejector  program  is  a  modification  of  Gilbert  and  Hill's 
ejector  analysis  program  (16).  It  includes  a  new  feature  which  can  design  a 
diffuser  with  the  incipient  separation  criterion  in  which  the  skin  friction 
velocity  is  maintained  at  a  low  non -zero  value  in  the  diffuser.  This 
criterion  produces  faster  pressure  recovery  in  the  diffuser,  and  therefore 
reduces  the  overall  ejector  length.  The  following  boundary  layer  type 
momentum  and  energy  equations  were  used  in  the  DTNSRDC  ejector  program  (17). 

X -Momentum 

3u  ur  3  I".  .  pur  3ul  1  3P 

"  jj*  2?3?L(,‘  +  £)  -Jruj-piJ  (  ) 

iji -Momentum 


Energy 


(3) 


(4) 


The  symbols  in  the  equations  are  defined  in  the  Nomenclature.  These  equations 
were  solved  using  the  finite  difference  method.  A  flow  diagram  for  the 
computer  program  and  instructions  for  input  of  data  are  reported  in  Appendix 
A. 

The  utilization  of  this  program  for  the  present  study  is  limited  to  the 
determination  of  the  velocity  and  temperature  distribution  up  to  the  exit  of 
the  mixing  chamber. 
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2.3  AVERAGE  VELOCITY  AND  TEMPERATURE  AT  THE  EXIT  OF  THE  MIXING  CHAMBER 


The  velocity  and  temperature  are  computed  using  the  DTNSRDC  Ejector 
Program  and  are  printed  out  for  all  the  selected  streamlines  at  preselected 
stations.  The  pressure  variations,  however,  are  only  considered  along 
streamlines  instead  of  across  streamlines.  Therefore  the  averaging  process 
only  deals  with  the  mass  averaged  velocity  and  temperature.  No  pressure 
averaging  is  necessary.  These  averages  are  given  by 


U) 


ave 


\updA 

=  TU) 


ave 


(5) 


and 


(T) 


S' 


TupdA 


ave 


A(p )  (u) 

ave  ave 


(6) 


2.4  DIFFUSER  AREA  RATIO,  CURVED-WALL  DIFFUSER  GEOMETRY  AND  THE  ROTATIONAL 
FLOW  CALCULATIONS 


2.4.1  Diffuser  Area  Ratio  —  Using  equations  (5)  and  (6),  the  Mach 
number  at  the  diffuser  inlet,  M^,  can  be  determined.  The  diffuser  inlet 
pressure  P^  is  merely  the  pressure  at  the  last  station  of  the  mixing 
chamber.  Assuming  the  flow  in  the  diffuser  is  isentropic,  and  the  exit 
pressure  Pe  is  the  ambient  pressure  for  the  analysis,  the  following  formula 
is  used  to  calculate  the  diffuser  area  ratio  AR. 


(7) 
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2,4.2  The  Curved-Wall  Diffuser  Geometry  —  The  geometry  of  the 
curved-wall  diffuser  is  determined  by  using  the  "Clemson  Inverse  Design 
Program  for  Short  Curved-Wall  Diffusers."  Reference  13  provides  the  necessary 
instructions  for  input  data  preparations.  For  plane-flow  cases,  Griffith 
diffusers  can  be  designed  by  specifying  only  the  area  ratio,  length-to-exit 
diameter  ratio,  and  the  percentage  of  flow  removal  used  for  boundary  layer 
control  (18).  For  axisymmetric  diffusers  the  current  version  of  the  computer 
program  requires  input  of  the  stream  function  and  the  velocity  potential  at 
the  boundary.  An  educated  guess  of  input  values  is  helpful. 


2,4.3  Rotational  Flow  Calculations  —  First  we  examine  the  complete 
vorticity  equation  in  vectorial  form  as  follows: 


(3(0  .“b  ->■  ->■  2-*- 

— - (-  (u  •  7)w  =  (u)  •  V)u  +  vV  iii 

<7  t 


(8) 


Consistent  with  the  approximations  made  in  this  study,  the  viscosity  term  is 
negle’  ‘ed  and  therefore: 


.-V  •+  — ►  -V  -K  -> 

—  +  (u  •  V)a,  =  (to  •  V)u 


(9) 


For  axisymmetric  flow,  the  term  (u)  •  V)  u  is  equal  to  — .  For  —  «1,  (^•'7  )u- 0 


and  consequently 


(10) 


Dm 

Dt 


or  vorticity  is  approximately  constant  throughout  the  diffuser. 


The  continuity  equation  for  incompressible  axisymmetric  flow  is: 


3v  3u  v 
3  r  3x  r 


0 


(ID 


Figure  6  defines  the  coordinate  system  and  the  velocity  components  used  in  the 
analysis. 
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The  stream  function  tj>(r,x)  is  defined  such  that: 


1  Hl 

u  =  —  r-1-  and  v 
r  3  r 


I  11 

r  3x 


and  then  equation  (11)  becomes: 


(12) 


2  2 
1  3ij)  _1  3  ij)  +  _1  3 

^2  3x  r  3r3x  r  3x3r 


1/1  3jA  = 

r  (  r  3x/ 


(13) 


In  addition,  ip  satisfied  the  vorticity  equation  which  in  cylindrical 
coordinates  for  an  axisymmetric  shear  flow  is: 


3v  3u 
3x  3  r 


(14) 


Substituting  equation  (12)  in  equation  (14)  yields 


3  1  |  | 

— o  “  ~  i  ~  w  r 

3 r2  r  3r 


(15) 


Which  is  an  elliptic  partial  differential  equation.  In  the  numerical  solution 
of  equation  (15)  for  rotational  flow,  we  use  the  same  grid  network  as  for 
irrotational  flow.  The  boundary  conditions  used  are  as  follows. 


2. 4. 3.1  Inlet  —  At  the  diffuser  inlet,  there  is  a  parallel  flow  in 
the  x-direction,  and  consequently  v  ■  0.  The  ip’s  are  obtained  from  equation 
(12)  as  follows 


1  3|  = 
r  3  r 


and 


_  Ill  =  0 

r  3x 


Therefore  <p  only  changes  with  r,  hence 


For  shear  flow, 


u 


ur 


u 

o ,  in 
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therefore 


Which  yields 


r  (u 


o,in 


wr)dr 


1  2 

tt  r  u  . 

2  o ,  in 


From  equation  (16),  we  have: 


wall 


iEi„ 


u 

o ,  m 


where  R  -  r(vWall> 


(16) 


(17) 


For  irrotational  flow  and  uniform  parallel  inlet  velocity  q^  ,  we  have: 


V  . 

irr 


/ 


in 


rq .  dr 
sirr 


At  the  wall: 


(v 


irr^wall 


1 

2  qirr 


(18) 


However,  may  not  necessarily  be  equal  to  (<{j  )  ,  which  means  that 

wall  irr  wall 

for  the  same  diffuser  geometry,  flow  rates  for  the  irrotational  and  rotational 

(shear)  cases  may  not  be  the  same. 


Define  Y  such  that 

^irr^wall  Yf,wall 


and  then 


q .  =  y  (u 

^irr  o 


.  -  4  ooR  ) 

,m  3  in 


(19) 
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For  the  shear  flow  solution,  we  will  need  to  (i)  specify,  for  the  inlet 
boundary  condition,  the  values  of  u  Q,in  anc^  (ii)  obtain  w  from  equation 
(19);  and  (iii)  specify  the  streamline  according  to  equation  (16). 

Another  possibility  is  that  both  w  and  u  0  are  specified,  and  Y  is  to  be 
calculated  from  equation  (19). 

In  the  inverse  design  program  if  we  have  prescribed  parallel  but  non-uniform 
irrotational  flows  at  the  inlet  and  the  exit,  then  the  boundary  condition 
at  the  inlet  would  be  as  follows: 

(it,  )  ,  -  yfr  R2  u  .  -  \  uR3  ] 

virr  wall  12  in  o,in  3  inJ 


For  known  u>  and  uq  y  can  be  calculated  from  the  following 

=  ^ir^wall 

{  R2  u  .  -7  wR3 

2  m  o,in  -  in 

and  the  streamlines  can  be  calculated  from  equation  (16). 

2. 4. 3. 2  Centerline 


(20) 


Along  the  centerline,  we  have  4*  *  0. 

2. 4. 3. 3  Wall,  upstream  of  slot 

Along  the  wall,  upstream  of  the  slot,  we  have  4*  *  ^wall  ’  w^ere  ^wall 
the  same  as  the  one  calculated  at  the  inlet. 

2. 4. 3. 4  Wall,  downstream  of  slot  and  inside  of  slot 

Let  B  and  Sirr  be  the  fractions  of  the  flow  into  the  slot  for  the 
rotational  and  irrotatonal  cases  respectively.  These  6's  are  related  to  the 
stream  functions  as  shown  in  the  following  expressions. 
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1 


<P 


eirr  1  ip. 


irr,st 


(21) 


i rr .wall 


and 


^wall  ^st 


wall 


or 


*st  ■  ‘wall'1  -  6) 


(22) 


3  is  specified  as  an  input  to  the  analysis.  The  ^  values  for  lines  AB  and  AC 
in  Figure  7  are  determined  from  equation  (21). 


2. 4. 3. 5  Diffuser  Exit 


Letting  Qirr  =  Volumetric  flow  rate  at  the  diffuser  exit  for 
irrotational  flow  and 


Q  =  Volumetric  flow  rate  at  diffuser  exit  for  shear  flow  results  in: 


Qirr/(1  ~  6irr}  =  yQ/(1  "  S) 


(23) 


It  can  be  shown  that  for  parallel  flow 


o.e 


T^stagu  + 


2(<p.  )  .,(1 

irr  wall 


-  8) 


yRNSTAGU 


(24) 


Here  w  is  specified,  '('irr  and  R^sTAGU  are  obtained  from  the  Clemson 
Inverse  Solution  Computer  Program,  and  uQ>e  is  the  center  line  velocity 
at  the  diffuser  exit.  (RNSTAGU  is  the  rac*ius  at  the  stagnation  point  A  of 
Figure  7.) 


2. 4. 3. 6  Exit  of  Slot 


At  the  slot  exit,  a  shear  flow  having  vorticity  w  is  assumed. 
Referring  to  Figure  8: 


R  -  R  =  y  sin(90  -  a)  =  y  cos  a 
o 
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J 


or 


R  =  R  +  y  cos  a 
o 


(25) 


The  following  development  is  to  determine  the  stream  function  in  the  slot. 

*y  ‘  *8t  =  JQ  qRdy 


(26) 


But  from  Figure  9  we  see  that,  q  «  9o,sl  “  ^ y 
and  therefore 

1  2  13 

to  =  to  +Rq  ,  y  +  —  (q  ,  cos  a  -  R  u>)y  -  -r  wy  cos  a 
y  st  o  o,sl  2  o,sl  o  3 


(27) 


The  'P's  at  the  exit  of  the  slot  can  be  prescribed  using  equation  (27)  where  y 

(,R-R  ) 

o 

is  replaced  by  cos  a  ,  vi2. 


(28) 


The  details  of  computational  methods  suitable  for  this  equation  subject  to  the 
proper  boundary  conditions  are  reported  in  references  19  and  20.  The  change 
of  flow  velocity  in  the  diffuser  is  much  less  than  that  in  the  mixing  chamber, 
therefore  in  this  part  of  the  computation  incompressible  flow  is  assumed,  even 
though  a  compressible  flow  approach  is  used  in  computing  the  diffuser  area 
ratio. 
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2.5  CRITICAL  VELOCITY  AND  SUCTION  FLOW  RATE 


Using  rotational  flow  analysis,  the  flow  velocities  U^  and  U2  outside 
the  boundary  layer  immediately  upstream  and  downstream  from  the  suction  slot 
are  determined.  From  Taylor's  criterion  (15)  the  critical  velocity  is 
determined  by  the  following  formula: 


The  following  expression  is  a  good  approximation  for  ucr 


u 


cr 


diff 


) 


2 


(30) 


where  AR^jj  is  the  diffuser  area  ratio. 


In  this  analysis  it  is  necessary  to  use  the  boundary  layer  velocity 
profile  at  the  exit  of  the  mixing  chamber  to  approximate  the  boundary  layer 
velocity  profile  immediately  upstream  of  the  suction  slot.  Therefore  the 
suction  flow  rate  can  be  determined  by  computing  the  mass  flow  rate  within  the 
boundary  layer  up  to  the  location  where  the  flow  velocity  is  equal  to  the 
critical  velocity  ucr.  This  suction  flow  rate  is  computed  as  follows 


R  -y 
me 


m  =  2ir 

suction 


purdr 


R 

e  ,mc 


(31) 


where  y  is  the  location  where  the  flow  velocity  is  equal  to  ucr. 
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2.6  GLOBAL  ANALYSIS  FOR  $2 

The  thrust  augmentation  ratio  $  2,  which  accounts  for  the  primary  flow 
in  the  auxiliary  ejector  and  the  thrust  contribution  of  its  discharge,  is 
defined  by 


m2V2  +  m3v3 


mV  +  m  , v  . 
P  P  P  P 


(32) 


where  m  is  mass  flow  rate.  The  subscripts  2  and  3  are  for  the  ejector  exit 
and  auxiliary  ejector  exit,  respectively,  subscript  p  is  for  the  primary  flow 
of  the  nozzle,  and  subscript  p'  represents  the  primary  flow  at  the  auxiliary 
nozzles. 


The  Global  analysis  was  performed  at  the  beginning  of  this  project  and  a 
one-dimensional  compressible  air-to-air  system  was  assumed.  The  primary 
purpose  of  the  Global  Analysis  was  to  examine  the  necessary  performance  levels 
of  mass  ratios  in  the  main  and  auxiliary  ejectors  in  order 
to  pre-select  a  value  of  j>2.  Toward  the  end  of  the  analytical 
investigation,  we  found  that  the  Global  analysis  is  useful  in  providing  a  ^ 
value  with  a  geometry  obtained  by  following  analytical  procedure  steps  1 
through  5  based  on  mass  averaged  velocity.  The  Global  analysis  is  attached  as 
Appendix  B. 

2.7  CORRECTION  OF  <f>2  FOR  NON-UNIFORM  VELOCITY 

Both  inlet  and  exit  of  the  diffuser  have  uni -directional  flows.  With 
known  vorticity  at  the  diffuser  inlet  and  the  constraint  of  constant  vorticity 
throughout  the  diffuser,  the  exit  flow  velocity  across  the  diffuser  can  be 
readily  computed.  The  relationship  between  a  mass-averaged  velocity  and  a 
momentum-averaged  velocity  can  be  calculated  for  a  known  non-uniform  velocity 
distribution,  thereby  providing  a  correction  for  $  2.  The  following  is  the 
formula  for  this  correction  factor, 
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3.  ANALYTICAL  PROCEDURES  FOR  INVESTIGATING  AN  EJECTOR  WITH  SHORT  CURVED-WALL 
DIFFUSER 

Figure  5  shows  the  major  steps  in  analyzing  an  axisymmetrical  ejector 
with  a  short,  curved-wall  diffuser.  The  mixing  chamber  is  shaped  like  a 
circular  pipe  with  a  contraction  toward  the  end  of  the  mixing  chamber.  To 
initiate  the  analysis  investigation,  the  area  ratio  of  the  primary  nozzle  to 
the  mixing  chamber  inlet  and  the  ratio  of  the  mixing  chamber  length  to  its 
diameter  (L/D)  are  specified.  The  diffuser  length  usually  is  restricted  to 
about  one  exit  diameter  of  the  diffuser.  In  addition,  the  static  pressure  at 
the  mixing  chamber  inlet  and  stagnation  pressure  of  the  primary  air  are  also 
specified.  Both  the  entrained  secondary  flow  from  the  stagnant  ambient  and 
the  primary  flow  within  the  nozzle  are  isentropic  up  to  the  mixing  chamber 
inlet.  Specifying  the  static  pressure  at  the  mixing  chamber  inlet  in  fact 
implies  specifying  the  ratio  of  the  secondary  mass  flow  to  the  primary  mass 
flow.  Starting  from  the  inlet,  there  is  heat  and  momentum  transfer  between 
the  primary  flow  and  the  entrained  secondary  flow. 

The  first  step  is  to  specify  the  mixing  chamber  geometry  by  a  selected 
value  of  C  along  with  other  mixing  chamber  geometric  parameters.  In  step  2, 
the  DTNSRDC  ejector  program  (17)  is  used  to  compute  velocity  and  temperature 
profiles  and  pressure  values  up  to  the  exit  of  the  mixing  chamber.  Governing 
equations  for  the  flow  are  of  the  boundary  layer  type;  therefore  pressure 


15 


variations  only  exist  along  the  flow  directions  and  not  across  the 
streamlines. 

In  the  third  step,  the  mass-averaged  flow  velocity  and  temperature  are 
computed.  In  the  fourth  step,  one-dimensional,  compressible  isentropic  flow 
is  assumed  in  determining  the  necessary  area  ratio  to  yield  the  diffuser  exit 
pressure  at  the  atmospheric  level  (or  the  ambient  level).  For  other 
computations  in  this  step,  inviscid  imcompressible  axisymmetric  flow  is 
assumed  within  the  diffuser.  As  long  as  the  flow  along  the  solid  wall  is 
maintained  without  deceleration,  the  inviscid  flow  approximation  can  be 
justified.  Under  the  framework  of  this  approximation  the  vorticity  of  the 
flow  will  remain  constant  everywhere  in  the  diffuser.  The  Clemson  Inverse 
Design  Program  (13)  is  used  to  obtain  the  geometry  of  the  short  curved-wall 
diffuser  for  a  specified  area  ratio  at  a  length  approximately  equal  to  one 
exit  diameter.  This  length  specification  is  somewhat  arbitrary,  yet 
experience  suggests  that  it  should  not  be  difficult  to  achieve.  Immediately 
after  the  determination  of  the  diffuser  geometry,  a  computer  program  from  the 
analysis  of  "Method  of  Rotational  Flow  Calculation  in  Griffith  Diffuser"  (19) 
is  used  to  examine  the  velocity  distribution  along  the  diffuser  wall.  It  is 
essential  that  there  be  no  deceleration  along  the  solid  diffuser  wall.  When 
necessary,  the  input  for  the  Inverse  Design  Program  is  revised  to  generate  a 
new  diffuser  geometry.  This  process  is  repeated  until  there  is  no 
deceleration  a^ng  the  diffuser  wall.  It  is  sometimes  necessary  to  revise  the 
analysis  from  the  very  first  step  where  the  mixing  chamber  inlet  pressure  is 
specified.  An  increase  of  the  static  pressure  at  the  inlet  implies  a 
reduction  in  mass  ratio  and  therefore  a  reduction  in  diffuser  area  ratio. 
Usually  this  revision  can  eliminate  the  problem  of  deceleration  along  the 
solid  wall.  For  a  configuration  with  no  deceleration,  no  flow  separation  will 
take  place.  It  is  recommended  that  the  velocity  from  the  diffuser  inlet  to 
the  suction  slot  be  kept  slightly  accelerated. 

In  order  to  estimate  the  amount  of  fluid  to  be  removed  by  boundary  layer 
control,  it  is  necessary  to  know  the  boundary  layer  profile  immediately 
upstream  of  the  suction  slot.  It  is  reasonable  to  assume  that  this  profile  is 
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the  same  as  that  at  the  exit  of  the  mixing  chamber,  and  by  using  that  profile 
to  calculate  the  critical  velocity  determined  by  Taylor's  criterion,  the  rate 
of  boundary  layer  removal  can  be  determined.  This  approximation  is  compatible 
with  other  approximations  used  throughout  this  analytical  procedure.  A  more 
accurate  computation  of  viscous  flow  within  the  diffuser  would  require  the  use 
of  the  Navier-Stokes  equation  because  appreciable  pressure  gradients  exist 
across  the  streamlines. 

It  is  apparent  that  the  quality  of  the  boundary  layer  at  the  mixing  cham¬ 
ber  exit  would  be  improved  if  the  mixing  chamber  had  a  contraction  toward  its 
exit  to  provide  acceleration.  With  a  contraction  a  smaller  amount  of  the 
fluid  within  the  boundary  layer  would  have  velocity  below  the  critical 
velocity.  On  the  other  hand,  a  contraction  of  the  mixing  chamber  will 
increase  the  diffuser  area  ratio  and  therefore  elevate  the  value  of  the 
critical  velocity  determined  by  Taylor's  criterion.  A  particular  combination 
of  these  two  counteracting  factors  will  lead  to  a  minimum  suction  requirement 
to  achieve  the  necessary  boundary  layer  control.  Therefore,  this  analytical 
procedure  is  formulated  with  mixing  chamber  contraction  as  the  optimization 
parameter. 

In  step  6,  the  modified  thrust  augmentation  ratio  4>  2  is  computed.  The 
consideration  of  the  mass  flow  of  primary  fluid  used  in  auxiliary  ejectors  to 
provide  the  necessary  boundary  layer  control  and  the  thrust  contribution  from 
the  discharge  of  the  boundary  layer  removal  are  included  in  the  definition  of 
4>  2  •  A  one-dimensional  compressible  flow  analysis  which  assumes  the  mass 
ratios  of  the  primary  and  auxiliary  ejectors  and  computes  the  mass  averaged 
velocity  at  various  parts  of  the  ejector  and  the  momentum  terms  at  the  ejector 
exit,  yields  the  thrust  augmentation  ratio  <f>2* 

In  step  3  the  velocity  across  the  mixing  chamber  exit  has  been  computed, 
and  in  step  4  an  inviscid  shear  flow  with  constant  vorticity  within  the  dif¬ 
fuser  has  been  approximated.  Therefore,  the  velocity  across  the  diffuser  exit 
can  be  readily  determined.  The  last  step.  Step  7,  is  to  correct  the  $2 
value  from  the  Global  Analysis  by  using  the  known  velocity  distribution  across 
the  diffuser  exit. 
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4.  RESULTS  AND  DISCUSSIONS 


4.1  COMPUTATIONAL  RESULTS 

The  following  computations  are  in  accordance  with  the  work  statements  of 
the  research  contract.  These  calculations  were  carried  out  for  the  purpose  of 
examining  the  validity  of  a  set  of  steam-to-air  ejector  data  and  to  provide 
specific  information  for  $  2  values  over  a  range  of  ejector  parameters  of 
immediate  interest. 

4.1.1  Steam-to-Air  Ejector  Analysis  Program  —  In  the  introduction 
section  of  this  report,  a  maximum  thrust  augmentation  ratio  of  2.81  was  cited 
as  an  estimated  value  from  the  experimental  mass  flow  ratio  and  velocity 
profile  data  of  reference  4.  The  analytical  tools  discussed  in  Section  2  were 
used  to  examine  the  experimental  data  of  reference  4  in  the  following  manner. 

4. 1.1.1  Computer  Program  Modif ications  —  It  is  recognized  that  the 
Hill  and  Gilbert  ejector  computer  program  (16)  was  developed  for  ejectors 
using  air  as  both  primary  and  secondary  fluids.  Modifications  in  computer 
programs  were  necessary  In  order  to  allow  steam  to  be  specified  as  the  primary 
fluid.  The  gas  constants  such  as  Y=  1.4  and  R  =  53.35  ft-lbf/lbm°R  for  air 
were  changed  to  y  =1.135  and  R=85.76  f t-lbf /lbm°R,  and  the  Prandtl  number 
was  changed  from  0.7  to  1.05  wherever  they  appear  in  the  computer  program  for 
calculations  involving  primary  flow.  The  turbulent  Prandtl  number  was  left 
unchanged.  In  addition,  in  the  modified  computer  program  the  identities  of 
streamlines  were  maintained  and  provisions  were  made  to  use  the  correct 
properties  (such  as  Pr)  for  each  streamline  in  the  finite  difference 
solutions. 


4. 1.1. 2  Comparison  of  Thrust  Augmentation  —  The  ejector  geometries 
used  in  this  comparison  have  mixing  chamber  inlet  diameters  and  nozzle 
discharge  diameters  the  same  as  those  used  in  reference  4.  Those  experiments 
were  conducted  with  nozzle  stagnation  pressures  of  10,  15,  20  and  30  psig,  and 
saturated  steam  was  the  primary  fluid. 
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The  computational  procedure  of  this  comparison  is  outlined  as 


follows: 


First,  the  inlet  radius  and  the  length  of  the  mixing  chamber,  the 
primary  nozzle  flow  area,  the  stagnation  pressure  and  temperature  of  the  pri¬ 
mary  and  secondary  air,  and  the  static  pressure  at  the  mixing  chamber  inlet 
are  specified.  Specifying  the  static  pressure  at  the  mixing  chamber  inlet  is 
equivalent  to  specifying  the  secondary  flow  rate.  The  values  of  static 
pressure  used  in  computations  were  taken  from  experimental  data  of  reference 
4.  After  these  values  were  specified,  the  DTNSRDC  ejector  design  computer 
program  was  used  to  provide  the  diffuser  geometry.  The  diffuser  obtained  from 
DTNSRDC  program  was  used  instead  of  Clemson  curved-wall  diffuser  because  the 
viscous  flow  is  assumed  to  be  of  boundary  layer  type.  The  geometries  of  Clem¬ 
son  diffusers  have  curvatures  too  large  to  be  considered  acceptable  for  a 
boundary  layer  flow.  A  computer  program  is  available  at  Clemson  to  accept  a 
curved-wall  diffuser  with  a  pressure  jump;  the  limit  on  the  pressure  jump, 
however,  is  too  severe  to  be  of  practical  value. 

Second,  the  modified  Hill  and  Gilbert  program  is  used  with  steam  as 
the  primary  fluid  and  air  as  the  secondary  fluid.  The  geometry  of  this  ejec¬ 
tor  is  identical  to  the  resulting  geometry  of  the  DTNSRDC  design  program.  The 
stagnation  temperature  and  pressure  of  the  primary  flow  have  the  same  values 
as  those  used  in  the  DTNSRDC  program.  The  static  pressure  of  the  secondary 
flow  at  the  inlet  of  the  mixing  chamber,  however,  is  no  longer  specified.  A 
new  mass  flow  ratio  is  obtained  as  a  part  of  the  output.  In  general  this  new 
mass  ratio,  using  steam  as  the  primary  flow,  is  about  30  to  40  percent  higher 
than  the  mass  ratio  of  an  air-to-air  ejector  obtained  with  the  DTNSRDC  design 
computer  program.  On  the  average  this  new  mass  ratio  is  about  20%  higher  than 
experimentally  obtained  values. 

Finally,  the  new  mass  flow  ratio  obtained  from  the  Hill  and  Gilbert 
program  is  used  as  input  to  the  global  analysis  program  to  compute  the  thrust 
augmentation  values  with  assumed  suction  rates  and  assumed  mass  ratios  for  the 
auxiliary  ejector.  In  the  global  analysis  the  thrust  augmentation  $2  is  de- 
fined  as  follows: 
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Figure  10  shows  three  sets  of  results  comparing  analyitcal  and  experimental 

values  of  thrust  augmentation.  Two  of  these  sets  represent  the  analytical 

results  of  thrust  augmentation  $2  versus  stagnation  pressure  of  the  primary 

flow  with  suction  flow  rates  of  5%  and  10%  and  the  mass  ratios  of  4,  5,  and  6- 

to-1  for  the  auxiliary  ejector.  The  remaining  single  curve  was  obtained  by 
•  •  • 

letting  m3  ■  mp  ■  0  and  m2  correspond  to  the  case  of  15%  suction  without 
a  secondary  ejector.  This  single  curve  can  now  be  compared  directly  with  the 
thrust  augmentation  values  <t>i  which  were  derived  from  the  measured  mass  flow 
ratio  and  the  exit  velocity  traverse  from  the  aforementioned  experiments  (4), 
viz . ,  using 


m2V2 

m  v 
P  P 


The  comparison  indicates  that  the  experimental  values  are  in  reasonably  good 
agreement  with  the  analytical  values. 


4.1.2 


Rotational  Flow  Calculations 


Figure  11  represents  an 


axisymmetric  diffuser  having  an  area  ratio  of  2.5-to-l.  Figure  12  shows  the 
velocty  distributions  along  the  diffuser  walls.  The  circles  denote  the 
velocity  distribution  prescribed  to  the  "Inverse  Design  Program”  for 
irrotational  flow,  and  the  triangles  represent  the  velocity  distribution  along 
the  diffusers  depicted  in  Figure  11.  The  second  velocity  distributions  were 
computed  for  inlet  shear  flow  having  a  nondimensional  vorticity  /alue  B  of 
0.33,  where  velocity  and  length  are  normalized  by  inlet  center  line  velocity 
and  inlet  width  respectively.  It  is  apparent  that  even  though  the  prescribed 
diffuser  wall  velocity  distribution  upstream  of  the  slot  has  a  strong 
acceleration,  the  wall  velocity  distributions  may  experience  acceleration  in 
the  diffuser  when  a  shear  flow  is  admitted.  In  actual  operation  a  flow 
separation  would  most  likely  take  place. 


Figure  13  shows  three  velocity  distributions  for  the  diffuser  depicted  in 
Figure  11;  these  are  specifically  for  inlet  non-dimensional  vorticity  values 
of  0.33,  0.5,  and  0.65.  There  is  a  significant  change  in  the  magnitude  of  the 
velocities,  but  their  gradients  vary  only  slightly.  As  expected,  larger 
vorticity  results  in  more  deceleration.  The  larger  change  of  wall  velocity 


magnitude  results  from  the  assumption  in  our  analysis  that  the  diffuser  center 
line  velocity  is  kept  constant,  therefore  the  wall  velocity  decreases  more 
when  the  inlet  shear  flow  vorticity  increases  more.  Judging  from  the 
distributions  of  this  figure,  a  flow  separation  would  likely  take  place  in  all 
three  cases.  Figure  14  shows  the  wall  velocity  distributions  for  the 
diffuser  depicted  in  Figure  11,  under  the  conditions  of  (i)  admitting  an  inlet 
shear  flow  with  a  vorticity  of  0.5  and  (ii)  three  different  suction  rates.  A 
suction  rate  of  5.6%  was  the  design  value.  At  a  suction  rate  of  8.5%,  the 
velocity  gradient  upstream  of  the  slot  has  not  improved  significantly  in 
comparison  to  the  distribution  resulting  from  the  5.6%  suction  flow.  Perhaps 
none  of  these  distributions  will  yield  a  fully  attached  flow  throughout  the 
diffuser. 

Figure  15  shows  a  diffuser  designed  for  irrotational  flow  and  having  an 
area  ratio  of  2.5-to-l  with  its  suction  slot  slightly  modified.  Figure  16 
shows  the  rotational  velocity  distribution  along  the  walls  of  diffusers 
depicted  in  Figures  11  and  15  with  inlet  vorticity  of  0.5.  Apparently,  there 
is  no  significant  difference  in  wall  velocity  distribution  for  these  two  shear 
flow  inlet  cases.  Figure  17  shows  a  diffuser  designed  for  irrotational  flow 
and  having  an  area  ratio  of  2-to-l.  Figure  18  shows  wall  velocity 

distribution  for  the  diffuser  dipicted  in  Figure  17.  The  velocity 
distribution  represented  by  circles  is  for  irrotational  flow,  and  that 
represented  by  triangles  is  for  shear  flow  with  vorticity  of  0.5.  In  this 
case  hardly  any  deceleration  is  detected  for  the  rotational  flow,  therefore, 
one  should  expect  a  highly  effective  diffuser  even  when  a  shear  flow  inlet 
condition  is  imposed  at  the  curved  wall  diffuser  inlet.  Figures  19  and  20  are 
similar  to  Figures  17  and  18  and  are  for  an  area  ratio  of  1.5-to-l. 

4.1.3  Results  from  Parametric  Study  —  The  parametric  study  of  thrust 

augmentation  ratio  <}>  2  for  an  air-to-air  ejector  was  made  for  a  set  of 

selective  values  with  MR“12,  (MR)  .  =  4,  F.S.  “  0.10  and  X  =  40.  The  mixing 

chamber  has  no  contraction  and  its  radius  is  1.4  inches.  The  diffuser  area 

ratios  are  around  2.0,  which  are  the  required  values  to  provide  the  exit 

pressure  level  to  match  the  ambient.  The  reasons  for  selecting  these  values 

are  given  in  the  section  4.2  of  this  report.  The  remaining  parameters, 

(L/D)  , ,  and  Pni/Pi*  were  each  varied  from  2.5  to  5.0  in  increments 

overall  1 
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of  0.5.  The  diffuser  length  was  designed  for  one  mixing  chamber  diameter 
approximately,  therefore  the  mixing  chamber  length  varied  from  1.5  D  to  4.0  D. 
The  ejector  inlet  static  pressure  was  slightly  below  the  atmospheric  pressure, 
therefore  the  nozzle  stagnation  pressure  varied  from  35  to  73  psia, 
approximately.  The  procedure  used  to  calculate  $2  is  descirbed  in  section 
3,  the  analytical  procedure. 

The  results  of  this  parametric  study  are  presented  in  Figure  21  as 
i>2  versus  P0[/p£  with  (L/D)  overall  as  parameter.  Figure  21  shows 
that  $2  increases  with  an  increase  of  the  pressure  ratio  and  also  increases 
with  a  decrease  in  length-to-diameter  ratio,  when  the  velocity  at  the  nozzle 
exit  was  assumed  to  be  sonic.  The  shaded  area  indicates  the  region  where  the 
vorticity  at  the  diffuser  inlet  is  too  high  for  rotational  flow  to  maintain  a 
constant  or  a  favorable  pressure  gradient  along  the  curved-wall  diffuser.  For 
a  short  curved-wall  diffuser  to  have  attached  flow,  this  requirement  must  be 
fulfilled.  The  computed  $2  value  will  no  longer  be  meaningful  when  the  flow 
within  the  diffuser  is  separated.  References  19  and  20  show  that  the 
non-dimensional  vorticity  should  be  kept  below  0.8  for  (AR)diff“2. 
Therefore  only  13  of  the  36  cases  computed  in  this  investigation  should  be 
considered.  Figure  22  shows  the  corrected  <p  2  values,  taking  into  account 
the  difference  between  mass-averaged  velocity  and  momentum-averaged  velocity. 
Only  the  $2  values  corresponding  to  those  13  cases  for  which  the  vorticity 
is  less  than  0.8  are  presented. 

Allowing  the  primary  flow  at  the  nozzle  exit  to  reach  supersonic  speed, 
the  <t12  values  will  reduce  to  the  levels  as  depicted  in  Figure  23.  The  trend 
of  thrust  augmentation  versus  pressure  ratio  exhibited  in  Figure  23  is  in 
agreement  with  the  traditional  trend. 

4.2  SAMPLE  EJECTORS 

Using  the  procedure  outlined  in  Section  3,  two  sample  ejectors  were 
designed.  These  two  ejectors  may  be  candidates  for  future  experimental 
investigation.  Figure  24  shows  one  of  the  sample  ejectors.  First  the  primary 
nozzle  diameter  was  selected.  A  diameter  of  0.448  inches  was  chosen,  based  on 
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the  capacity  of  a  compressed  air  system  in  the  Mechanical  Engineering 
Department  at  Clemson  University. 

Next  the  design  specification  X  for  the  area  ratio  of  the  mixing  chamber 
to  the  primary  nozzle  was  selected.  This  ratio  is  somewhat  limited  by  the 
amount  of  space  available  in  a  particular  application.  Since  this  thrust 
augmentation  may  have  application  in  VSTOL  aircraft,  relatively  smaller 
values,  specifically  40  and  20,  were  used  in  comparison  with  the  X  values  of 
reference  4. 

Based  on  the  above  specifications,  the  mixing  chamber  inlet  diameters  are 
2.8  inches  and  2.0  inches  for  X  =  40  and  20,  respectively,  and  each  mixing 
chamber  has  a  10%  area  reduction  for  contraction  at  its  exit.  The  ratio  of 
length  of  the  mixing  chamber  to  its  inlet  diameter,  L/D,  was  selected  to  be  4. 
This  length  is  not  an  optimum  choice  but  is  judged  to  be  adequate  and  short 
enough  not  to  be  deemed  excessive  when  considering  viscous  dissipation.  A 
more  adequate  discussion  will  be  given  in  this  regard  in  the  section  4.3  on 
the  optimization  process. 

The  diffuser  area  ratios  are  2.187  and  1.993  for  X  *  40  and  20, 

respectively.  The  corresponding  ratios  of  overall  length  to  mixing  chamber 
inlet  diameter  are  5.256  and  5.182,  and  mass  ratios  of  12-to-l  and  7-to-l  are 
predicted  with  specified  mixing  chamber  inlet  static  pressure  levels. 

The  corresponding  $2  values  with  the  corrections  for  non-uniform 

diffuser  exit  velocities  are  predicted  to  be  1.50  and  1.30  for  X  =  40  and  20, 
respectively,  providing  that  a  mass  ratio  of  4-to-l  is  accomplished  in  the 
auxiliary  ejectors  for  boundary  layer  control.  In  the  literature,  A.  E.  Kroll 
(3)  reports  that  a  mass  flow  ratio  of  3  to  5-to-l  can  be  achieved  with 

conditions  similar  to  those  of  the  auxiliary  ejector.  Using  the  DTNSRDC 

Ejector  Program  with  the  above  mixing  chambers,  these  values  of  mass  ratio, 
12-to-l  and  7-to-l,  were  shown  to  be  possible  with  diffusers  designed  with  the 
incipient  separation  criterion. 
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4.3  OPTIMIZATION 


Using  the  mixing  chamber  contraction  to  optimize  the  thrust  augmentation 
ratio  was  discussed  in  Section  3  on  the  analytical  procedure.  Figure  25  shows 
the  fraction  of  suction  required  for  necessary  boundary  layer  control  for  MR  * 
12  and  \  »  40.  Figure  26  shows  the  estimated  thrust  augmentation  ratios  <f>2 
for  these  specifications.  The  <j>  2  increases  slightly  from  0%  contraction  to 
10%  contraction.  For  contractions  of  10%  to  15%,  the  <t>2  values  are 
virtually  constant,  but  a  rapid  decrease  occurs  beyond  20%  contraction.  A 
greater  degree  of  contraction  will  produce  a  larger  area  ratio  for  the 
diffuser  and  may  make  it  difficult  to  maintain  non-deceleration  along  the 
diffuser  wall.  Therefore  the  optimum  contraction  is  around  10%.  The  results 
for  ejectors  of  MR  =  7  and  20  are  similar.  A  slight  increase  in  $2  was 
observed  from  0%  to  10%  contraction.  Beyond  15%  contraction,  decreases  in 
<f>2are  observed  due  to  a  rapid  increase  in  the  required  suction  flow  rate. 

Optimization  also  can  be  applied  to  the  mixing  chamber  length  in  the  same 
way.  It  was  found  that  for  MR  *  12  and  1  =  40,  ^  would  increase  from  1.50 
to  1.63  if  the  mixing  chamber  length  is  8.4  inches  instead  of  the  current 
design  of  11.2.  The  reduction  of  the  mixing  chamber  length  will  not  only 
result  in  a  decrease  of  boundary  layer  thickness,  but  also  cause  an  increase 
of  vorticity  of  the  shear  flow  at  the  diffuser  inlet.  The  second  effect  makes 
the  maintenance  of  non-decelerating  wall  velocity  more  difficult.  These  two 
counteracting  factors  make  the  optimization  process  for  length  possible. 


5.  CONCLUSIONS 


The  major  conclusions  which  were  obtained  from  this  investigation  are: 

1.  Analysis  shows  that  a  significant  increase  in  the  thrust  augmentation 
ratio  of  the  thrust  ejector  can  be  achieved  with  significantly  reduced  overall 
length  when  a  Griffith  diffuser  is  used.  Auxiliary  jets  to  provide  the 
necessary  boundary  layer  control  seem  attractive  because  they  are  compact  and 
simple  in  structure.  Their  effectiveness  as  a  means  of  boundary  layer  control 
remains  to  be  demonstrated  experimentally. 


2.  In  contrast  to  the  irrotational  flow  case,  for  each  value  of  vorticity 
there  is  an  upper  limit  for  the  area  ratio  within  which  the  diffuser  can  be 
designed  without  wall  deceleration. 

3.  Reducing  the  diffuser  area  ratio  is  the  most  effective  way  to  avoid 
deceleration  along  the  diffuser  wall.  This  implies  that  the  initially 
specified  mass  flow  rate  of  the  primary  ejector  must  be  reduced. 

4.  Both  the  mixing  chamber  length  and  a  contraction  at  the  mixing  chamber 
exit  can  be  used  as  optimization  parameters  for  thrust  augmentation  ratio. 
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Figure 


1.  A  sketch  of  a  thrust  augmentor  with  short  curved- 
wail  diffuser.  Secondary  jets  are  used  to  provide 
the  necessary  boundary  layer  control. 


Typical  velocity  distribution  in  a  Griffith 
Diffuser. 
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Figure  3.  Typical  wall  contour  of  a  Griffith  Diffuser. 


Figure  4.  Typical  velocity  variation  along  the  streamlines 
across  the  diffuser 
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Figure  5 
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"DTNSRDC  Program" 
Computation  up  to 
mixing  chamber  exit 


Average  temperature  and 
velocity  and  pressure 
at  mixing  chamber  exit 

Compute  diffuser  area  ratio  AR. 
Go  to  "Inverse  Design  Program" 
and  "Shear  Flow  Calculation"  | 
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"Global  Analysis" _ 


Correct  <t>2  for 
non-uniform  flow 


Block-diagram  showing  the  major  steps  in  ejector  analysis 
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Figure  6.  Coordinate  system  and 
velocity  components 
used  in  analysis. 


Figure  7.  Wall  streamlines  and 
the  branching  point. 


Figure  8.  Coordinate  system  at  slot  exit. 


Figure  9.  Velocity  distributions  for  shear  flow. 
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Figure  21.  Thrust  augmentation  ratio  tf>2  vs.  oressure  ratio  wit! 
(L/D)overall  35  Par£meter- 


Thrust  augmentation  ratio  4>2,  corrected 


Figure  22.  Corrected  thrust  augmentation  ratio  412  vs.  pressure  ratio 
with  (L/°)overall  as  parameter. 
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Corrected  thrust  augmentation  ratio  i>_  vs  pressure 
ratio  with  supersonic  primary  flows  and  with 
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Figure  24.  Sample  ejector  (supersonic  V  and  V  , ). 
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Figure  25.  Boundary  layer  removal  vs.  Figure  26.  Thrust  augmentation  ratio  <f>2 

percentage  of  contraction.  vs.  percentage  of  contraction. 
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APPENDIX  A 


In  the  DTNSRDC  program,  ejector  inlet  conditions  are  prescribed,  as  well 
as  the  geometry  of  the  mixing  chamber  of  the  ejector.  The  inlet  conditions 
are  the  relative  dimensions  of  primary  nozzle  exit  and  mixing  section  inlet, 
and  the  relative  mass  flow  rates  of  the  primary  and  secondary  flows.  The 
prescribed  geometry  is  used  only  up  to  the  point  of  application  of  the 
incipient  separation  criterion,  corresponding  to  the  diffuser  inlet. 

Figures  A-l  and  A-2  represent  flow  charts  for  the  overall  program  and  the 
inner  loop  for  the  program.  The  inner  loop  contains  a  more  involved  iterative 
procedure  than  was  present  in  Gilbert  and  Hill's  program;  the  pressure  gradi¬ 
ent  in  the  diffuser  section  (REXO  >  0  or  X  _>  XO)  is  adjusted  iteratively  to 
meet  the  incipient  separation  criterion. 

The  program  contains  two  termination  tests: 

1)  If  the  number  of  iterations  in  the  inner  loop  gets  larger  than  99  without 
the  incipient  separation  criterion  being  satisfied,  or  if  the  friction 
velocity  is  negative  at  the  end  of  the  inner  loop,  abnormal  ending  of  the 
calculations  occur. 

2)  If  the  maximum  length  prescribed  is  reached  or  if  atmospheric  pressure  is 
achieved,  normal  ending  of  the  calculations  occurs. 

The  following  is  the  set  of  instructions  for  data  input. 
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INSTRUCTIONS 


for  data  input 


Card  No. 
1 
2 

3 

4 
3 
6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 


Parameters 

DP  10 ,  DP29,  DP11 

MR,  NS,  N,  NJ,  NK 


Format 
6E13.6 
815,  4F10.8 


PROFE(I),  I  "  l»  MK 


8F10.8 


RRO(I) ,  I  =  U  NS 


8F10.8 


XXO(I),  I  =  i.  NS 


8F10.8 


PSO(I) »  I  *  53’  70 


6E13.6 


NP,  NPP,  NL, 
IPUNCH,  KJ, 


NTEST,  NJ,  NM 
SQP,  CL,  CU  PERCNT 


POl,  TOl ,  PAMB,  TOO,  PI, 
YJOU,  RD 


815,  4F10.8 
8F10.8 


17 


X,  XX,  XQ,  REXO,  CP2,  CXX2, 
CA,  B 


8F10.8 


18 


Card  18  is  the 
must  be  set  a 


NP,  NPP,  NL, 
IPUNCH,  KJ, 


NTEST , 
SQP,  CL 


NJ,  NM, 

,  CU,  PERCNT 


8 L 5  ,  4F10.8 


last  card  to  end  the  calculations 
value  larger  than  1. 


of  the  program, 


on 


which  NPP 
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Cards  1  through  17  are  required  for  each  set  of  data.  For  data  more  than  one 
set,  cards  1  through  19  must  be  placed  in  the  same  sequence. 


DP10 ,  DP 20 ,  DP  1 1 


MK 

NS 

N 

NJ 

NK 

PROFE(I) 


BRO(I) 

XXO(I) 

PSO(I) 

NP,  NPP ,  ML,  SQP 

NTEST 
CL,  CU 


POl 


MEANINGS  OF  SYMBOLS 

Initial  guessed  pressure  gradients  on  m  *  1,  2  and 
zero  lines  respectively.  The  initial  guesses  of 
the  values  of  DP10,  DP20,  DP11  at  the  initial  plane 
may  be  assumed  equal  at  any  plus  or  minus  dimen¬ 
sionless  value  of  the  order  of  10”^  to  10“®. 

A  number  which  indicates  the  number  of  velocity  and 
temperature  detail  being  printed  out. 

Indicates  the  total  number  of  stations  at  which  the 
geometry  is  given  in  sub-routine  RADIUS. 

Total  number  of  streamlines. 

Number  of  streamlines  for  primary  flow. 

Number  of  streamlines  for  secondary  flow,  outside 
the  boundary  layer. 

An  array  containing  MK  values  of  downstream  posi¬ 
tions  in  inches  where  the  velocity  and  temperature 
details  are  required  to  be  printed  out. 

An  array  containing  the  total  number  (NS)  of  duct 
radius  (inches). 

An  array  containing  the  total  number  NS  of  axial 
stations  where  the  RR0(l)'s  are  provided  (inches). 

Last  18  values  of  the  non-dimensional  stream  func¬ 
tion  specified  to  satisfy  the  boundary  layer  re¬ 
quirements  (see  reference  1  for  sample). 

Control  variables.  For  axisymmetric  flow,  NP  *  2, 
NPP  =  0,  NL  -  1,  SQP  =*  0.5. 

Test  run  number. 

Used  in  applying  the  incipient  separation  cri¬ 
terion.  They  are  lower  and  upper  limits  to  the 
dimensionless  velocity  at  the  last  streamlines  be¬ 
fore  the  wall.  CL  ■  0.0035  and  CU  ■  0.0040  are 
representative  values. 

i 

Nozzle  plenum  chamber  pressure  (PSIA). 
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T01 

PAMB 

TOO 

PI 

YJIN 

YJOU 

RD 

X 

XX 

XO 

REXO 

CP2 


Nozzle  plenum  chamber  temperature  (°R). 

Ambient  pressure  (PSIA). 

Ambient  temperature  (°R). 

Pressure  at  nozzle  exit  (PSIA). 

Nozzle  exit  inner  radius . 

Nozzle  exit  outer  radius. 

Mixing  section  radius  at  inlet. 

Distance  from  duct  inlet  to  nozzle  exit  plane  at 
which  calculation  begins  (inches). 

Distance  from  duct  inlet  at  which  calculation 
stops . 

Distance  from  duct  inlet  at  which  incipient 
separation  criterion  is  applied. 

Control  variable.  Is  negative  for  X  >  XO,  positive 
for  X  >  XO.  Set  REXO  «  -1  for  input. 

Gauge  pressure  (dimensionless)  at  diffuser  exit. 
Set  CP 2  =  0. 


CXX2,  CA,  B  Parameters  used  in  pressure  gradient  calculation 

for  incipient  separation  criterion.  Set  CXX2  * 
1.005,  CA  =  0.325,  B  «  1. 
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Flow  Chart  of  DTNSRDC  Computer  Program 


CALCULATE  jj> 


SI 


ites  >  9 


IS  'OUT  3AT 1 SFACT3RY 


.‘10  YES  ! 


*E5ET  mEW  vAluEs  it  Y,  — 
3  equal  to  old  stored  ones 


Figure  A-2.  Inner  Loop  of  DTNSRDC  Computer  Program 
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APPENDIX  B 

GLOBAL  ANALYSIS  FOR  PARAMETRIC  STUDY  OF 
THRUST  AUGMENTING  EJECTOR  USING  GRIFFITH  DIFFUSER 
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1 


Introduction 


This  analysis  is  conducted  for  the  purpose  of  determining  the  ranges  of 
operating  parameters  which  will  lead  to  a  thrust  ratio  above  2.0  in  an 
ejector  for  thrust  augmentation.  The  thrust  ratio  is  defined  as  the  ratio 
of  the  exit  momentum  of  the  mixed  fluid  to  the  momentum  of  the  primary  fluid 
at  the  inlet. 


In  this  study  a  special  type  of  curved  well  diffuser,  referred  to  as  a 
Griffith  Diffuser,  is  used.  In  order  to  operate  the  diffuser  at  a  very  high 
efficiency  as  it  is  intended,  boundary  layer  control  is  necessary.  Several 
auxiliary  jets  will  be  used  to  achieve  the  required  boundary  layer  control. 
A  parametric  study  was  conducted  in  which  net  thrust  rations  were  computed 
in  addition  to  the  thrust  ratios  defined  above.  The  net  thrust  ratio  is 
defined  as  the  results  of  dividing  the  sum  of  exit  thrust  and  pressure 
thrust  minus  the  ram  drag  by  the  momentum  of  the  primary  flow.  Upon 
examining  the  net  thrust  ratio  one  may  see  the  relative  contributions  from 
the  exit  momentum  terms  and  the  pressure  terms  as  well  as  the  liability  of 
the  ram  drag  attributed  by  the  momentum  at  the  inlet  of  an  ejector.  Figure 
B1  shows  a  schematic  diagram  of  the  ejector  for  this  study. 


- 


Figure  BI.  Schematic  diagram  of  an  ejector  for  thrust 
augmentation.  Auxiliary  jets  are  used  for 
boundary  layer  control  in  diffuser. 

A  A 

(A,  “A  +  A  A  ,  «  A .  -  A  -  JsA  and  -p—  "  - — ^ — —  ) 

pi  p  p  si  1  p  p 


si 


A  -  A 
1  P 
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\  1 


—  -  » 
-  /  V-J 


;-q  *  \ 


K  A, 


Figure  B2.  Control  volume  for  momentum  survey. 


Analysis 


Figure  B2  shows  a  control  volume  from  which  the  momentum  equation  for 
the  fluid  in  x-direction  is  written  as  follows. 


P  A  +  p  A  -  p„  A  -  p  A  +  Th 
Pi  Pi  si  si  2  2  3  3  fx 

*  1  [p  V  +  A„  V  -  A  V  -  A  V  ] 
7-  2  2  33  pp  si  si 

sc 


where  Thfx  ■  force  component  in  x-direction  on  the  fluid  at  the 
Interface  of  ejector  and  the  fluid  inside  the  control  surface. 

Ih  •  1  !A  V„  +  A„  V  -  A  V  -  A  V  ] 

fx  —  '■  2  2  33  pp  »i  si 

8c 

-  P  A  -  p  A  +  p,  A  +  ? ,  A 
Pi  Pi  Si  sx  22  33 

The  sum  of  force  components  in  the  x-direction  on  the  ejector  is  denoted  by 
Th*jx  and  consists  of  the  forces  over  the  interior  surface  and 
exterior  surface.  The  former  is  equal  and  opposite  to  Thfx.  Therefore, 
Th€j  -  -  Thfx  +  force  components  in  x-direction  on  the  exterior 
surface. 


I 


To  find  the  force  components  in  x-direction  on  the  exterior  surface  of  the 
ejector  we  may  consider  the  following: 


Figure  B3.  Exterior  forces  over  the  ejection. 


Figure  B3  shows  the  pressure  force  system  over  the  exterior  surface  of  the 
ejector  and  can  be  replaced  by  an  equivalent  of  two  sub-systems.  Figure  B4 
shows  the  sub-systems. 


Figure  B4.  An  equivalent  force  system  made  up  of  two  sub-systems. 


The  left-hand  side  sub-system  has  a  resultant  of  zero,  therefore  the  force 
system  represents  the  exterior  pressure  force,  shown  as  Figure  B3,  can  be 
replaced  by  the  right-hand  side  of  the  sub-system  of  Figure  B4  or 


-p  A  +  p  A_  +  p  A,  +  p  A  +p  A 
^atm  ‘atm  3  atm  2  pi  pi  9i  si 


-  >2  A2  -  p3  *3  '  Jc  [42  72  +  43  ?3  *  4p  Vp  '  4s,  T.,] 

Since  P2  and  P3  are  the  ambient  pressure,  we  have  p^  *  P3  “  Patm* 
In  addition, 


Therefore, 


Th  .  -  -  p  „  (A  +  A  )  +  p  A  +p  A 

ejx  atm  sj.  Pi  *Pi  Pi  «i  sl 

-  -  6a-  V  +  a ,  V  -  m*  V  -  m  V  ] 


2  2 


3  3 


P  P 


si  si 


Since  the  primary  jet  is  placed  away  from  the  inlet  of  the  mixing  chamber, 
we  have  ppl  -  psl 


Therefore, 


Th 


ejx 


(P. 


-  p  )  A  -  —  (m  V  +  m  V.)  +  —  (m  V  +  m  V  ) 
atm  i  g  2  2  3  J  gc  P  P  si  si 


In  this  equation,  (Psl  “  patm)  is  th®  gage  pressure  at  the 

mixing  chamber  inlet  and  is  a  negative  quantity  and  contributes  positively 
to  the  thrust  on  the  ejector.  We  refer  to  this  as  pressure  thrust.  The 
term  [-  —  (*2  ^2  +  m3  V3)  ]  is  the  exit  thrust  force  that  also 
contributes  positively  to  the  thrust  on  the  ejector.  The  last  term, 
[(  i  (ip  Vp  +  ®si  Vs  j)]  is  the  ram  drag,  a  liability  to  the 

thrus'l  of  the  ejector.  If  the  primary  jet  is  located  inside  the  mixing 
chamber  inlet,  then  the  condition  ppj  *  pgj  r.«iy  not  hold. 


In  a  single  thrust  augmentor,  the  thrust  ratio  is  usually  denoted  by 
and  is  defined  as  follows: 


♦ 


m  V 
P  ? 


In  the  system  described  in  this  analysis  we  use  the  symbol  $2  f°r  thrust 
ratio,  and  this  is  defined  as  follows 


°2  V2  *  °3  V3 
m  V  +  a  ,  V  , 
P  P  P  P 


where  «p*  and  Vpt  denote  the  mean  flow  rate  and  jet  velocity, 
respectively,  of  the  auxilliary  jets  used  in  the  diffuser  for  boundary  layer 
control. 
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The  net  thrust  ratio  then  is  denoted  by  TR,  and 


Th 
TR 


(m  V  +  m  ,  V  , )  i 
?  P  P  P  g„ 


-On.  v  +  m  V  )  +  m  V  +  m  V  +  (p  )  Ag 
2  2  3  3 _ p  p  si  si  si  gage  i  c 

m  V  4-  m  V 
P  P  P  P 

For  convenience  we  introduce  the  following  symbols: 


MR  *  mass  ratio  *  t 


m  m  , 

3  :  and  MR  -  -r 2_ 


m 


p  m 


Then 


in  •  MR  a  and 

si  p 


m,  -  (1  -  FS) (MR  +  l)m 
*•  c 


where  FS  stands  for  fraction  of  through  flow  sucked  away  for  the  purpose  of 
boundary  layer  control  in  the  diffuser,  and 

m_  *  m  ,  +  FS  (MR  +  1)  m  . 

3  p  p 


It  can  be  shown  that  the  net  thrust  ratio  TR  can  be  expressed  as  follows. 


-(1  -  FS)  (MR  +■  1)V,  -  MR  V  -  FS  (MR  +■  1)V 
TR  ■  -  P  3  3 

V  +  MR  V  . 

P  P  P 

(exit  thrust  ratio) 
g„ 

(p  )  A  tS- 

si  gage  i  a 

+ _ l _ 

V  +  MR  V  ,  + 

P  P  P 


(pressure  thrust  ratio) 


V  +  MR  V 

_£ _ ll_ 

V  +  MR  V  , 
P  P  P 
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(ram  drag  ratio) 


Introducing  PR,  the  absolute  pressure  ratio  of  no22le  plenum  to  secon¬ 
dary  flow  at  the  mixing  chamber  inlet 


PR  = 


o,a 


si  ,a 


Then  (pgi) 


gage 


o  ,a 
PR 


-  P 


atm 


and 


TR 


V 

V 

P 


+  (MR)  V 

_ s_ 

+  (MR)  V  , 
P  P 


(111* 

K  PR 


-?atm)(1+r 


0  V 

JLJL 


(MR)  V 
P  P  i 


-  ^2 


Geometrical  and  Flow  Parameters 


The  ratio  of  the  secondary  flow  area  at  the  ejector  inlet  to  that  of 

,  and  the  diffuser  area  ratio,  (AR)diff,  are 

The  mixing  chamber  geometry,  however,  at  this 
stage  of  the  study  will  not  be  specified  directly.  For  each  mixing  chamber 
along  with  other  parts  of  the  ejector  under  a  particular  set  of  operating 
conditons  (namely  the  plenum  pressure  of  the  primary  jet  P0>a  and  the 
discharge  pressures  P£  and  P3),  there  is  a  specific  velocity  profile 
within  the  boundary  layer  upstream  of  the  suction  slot  of  the  diffuser. 
This  specific  boundary  layer  implies  a  specific  value  of  FS ,  the  fraction  of 
suction,  and  a  specific  value  of  mass  ratio  MR.  Therefore  we  will  assign  a 
range  of  values  for  FS  and  MR  respectively  to  suggest  that  different  mixing 
chambers  are  involved. 

Similarly  we  will  assign  several  values  of  MRp,  the  ratio  of  mass 
flow  rate  in  the  primary  jet  to  that  in  the  auxiliary  Jets,  instead  of  spe¬ 
cific  geometrical  details  of  the  auxiliary  jets.  Alternatively 
(MR)aux,  the  ratio  of  entrained  mass  flow  rate  to  the  mass  flow  rate 
in  the  jets  of  the  auxiliary  ejector,  for  boundary  layer  control  may  be 


the  primary  jet,  .  S1 

A 

P 


specified  geometrically. 
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introduced  and  MRp  calculated  from  (MR)aux  and  other  parameters. 
Specifying  (MR)aux  along  with  MR  provides  a  relative  performance 
description  between  the  primary  ejector  and  the  auxiliary  ejector. 

Once  the  above  parameters  are  specified,  all  the  variables  involved  in 
computing  <t>0  and  TR  can  be  determined  from  the  following  equations. 


1. 


V 


P: 


V 

P 


2g  J  1  C  T  [1 
c  N  p  o 


k-1 

*  ] 


1_ 

2 


where  jj  is  the  nozzle  efficiency,  and  is  assumed  to  be  0.95 
Cp  is  the  specific  heat  at  constant  pressure  for  air 
k  is  the  ratio  of  specific  heats 

p_  _  is  either  the  atmospheric  pressure  or  the  critical 

rPla 

pressure,  0.5283  P0>a 

P0  a  is  previously  defined  and  is  specified 
T0  is  air  temperature  inside  the  nozzle  plenum 


2. 


P‘ 


B  >a_ 


RT 


pa 


P  k-1 
R  T  -2^-  ~ 
°  po,a 


where  R  is  the  gas  constant 


3.  V  : 
si 

From  m  *  MR  m 

s  p 


(p  A  V)  -  MR  (p  A  V) 
si  P 


-  (MR)  (DR)  (\)  Vp 
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r 


where  DR 


a  cm  d  ,a 


acm 


A 

C 

a” 


o  ,a 


k- 1 
k 


(Ps) 


gage 


gage  -g. 


5.  7p.: 


-g,  J  n. 


C  T 
N  p  o 


Ll  - 


o  ,a 


For  the  range  of  P0>a  to  be  considered 


Pp  ”  Pp'  *  Pcritical 
Therefore,  Vp  »  Vpi 


6. 


MRp: 

MR? 


a  A 
P  P 


V 


V 

? 


7. 


(MR)aux: 


m  ,  (MR)  »  F.S.  (MR  +  l)m 
p  aux  p 


(MR) 


aux 


F.S.  (MR  +  1) 
MR 

P 
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8.  V2: 


(pAV),  *  (m  +  m  )(1  -  F.S.) 

2  si  p 

(MR  +  1) (1  -  F.S.)  P 

.  ’ .  V  * - E. 

2  (AR)  A!  02 


-  (MR  +  1) (i  -  F.S.  )m 

P 


A  V 

_P—  E. 


where  p2  is  at  T2  and  P2  corresponding  to  standard 
pheric  conditions. 


9.  V3: 

<®3  v3>x  *  <V  VPf)x 

where  the  subscript  x  denotes  x-direction 
and  m3  *  F.S.  (MR  +  1)  thp  +  Api 
(MR) 

.  *  V  =  - —2 -  V 

3  F.S.  (MR  +  1)  +  MR?  p’ 


Range  of  Values  Used  in  Parametric  Study 

1.  AD  =  0.0013  ft2 

P  A  l 

2.  X  *  -  103.9,  176.3 

A 

3.  AR  -  ?.0,  1.25,  1.50 

4.  F.S.  -  0.05,  0.10,  0.15 

5.  MR  -  10,  13,  16,  19,  22,  25 
6-  (MR)aux  _  5j  7.5 

7.  (PQ)  -  35,  43.75,  52.5,  61.25,  70  (PSIA) 


atmos- 
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Results:  One  thousand  and  eighty  (1080)  TR  and  <$i  were  computed. 

68%  of  the  1080  runs  led  to  2  above  2.0.  The  results  of  this  analysis 
ponts  out  the  necessary  levels  of  performances  for  MR  and  (MR)aux  and 
th  restrictions  on  boundary  layer  development  in  the  mixing  chamber  required 
to  achieve  a  thrust  ratio  above  2.0. 


